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Abstract 

It is well known tinat the Hankel transform possesses neitlier a shift-modulation nor a convolution-multiplication 
rule, both of which have found many uses when used with other integral transforms. In this paper, the generalized 
shift operator, as defined by Levitan, is applied to the Hankel transform. It is shown that under this generalized 
definition of shift, both convolution and shift theorems now apply to the Hankel transform. The operation of a 
generalized shift is compared to that of a simple shift via example. 



Introduction 

It is well known that the Hankel transform does not satisfy 
a convolution or shift theorem in the simple way as that 
the Fourier and Laplace transforms (Piessens 2000), re- 
ducing its apparent utility. This follows because the 
Bessel functions do not possess a simple addition formula 
in the same way that the exponentials satisfy e'*^^-''' = e^'e'-'. 
The operation of convolution between two functions 
consists of a shift in one of the functions, a multiplica- 
tion with the other function, followed by an integration 
(or summation for discrete transforms) over all allow- 
able shifts. Thus, the lack of a convolution theorem for 
the Hankel transform follows because of the lack of a 
simple expression for the shift of a function in the 
Hankel transform domain. 

Levitan introduced the idea of a generalized displace- 
ment operator (Levitan 2002). As useful as this concept 
might be, to the best of the author's knowledge it does not 
appear to have seen much application in the physics or en- 
gineering literature. In this paper, Levitan's generalized dis- 
placement operator is made concrete via application to 
the Hankel transform. We show that this leads to shift 
and convolution rules for the Hankel transform. By way of 
several examples, the Hankel (generalized) shift is com- 
pared to the standard simple shift. 



Background 

Several definitions of the Hankel transform appear in the 
literature. In this paper, we use the definition of the «th 
order Hankel transform as defined by Piessens in (Piessens 
2000) to define the Hankel transform as 



F{p) = H„(/(r)) = J f{r)h{pr)rdr, (1) 

0 

where /„(z) is the «th order Bessel function. Here, n may 
be an arbitrary real or complex number. However, an inte- 
gral transform needs to be invertible in order to be useful 
and this restricts the allowable values of «. If « is real and 
n> - 1/2, and under suitable conditions of integrability of 
the function, the transform is self-reciprocating and the in- 
version formula is given by 

/(r) = J F{p)J„{pr)pdp (2) 

0 

Thus, the Hankel transforms takes a function j{r) in 
the spatial r domain and transforms it to a function F(p) 
in the frequency p domain. This relationship is denoted 
symbolically asfir)<>F{p). 
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Generalized shift for the IHankei Transform 

(Levitan 2002) introduced the idea of a generalized shift 
operator. Using Levitan's definition, the generalized shift 
operator R"^" indicates a shift of ro acting on the function 
J{r) is defined by the formula (Levitan 2002). 



where the order of integration has been reversed 
and where: 



(6) 



D„ (r, ro,x) = / /„ (^r)/„ {pro)Jn {px)pdp 



R^ir) =/(r|ro) = / F{p)J„{pr,)J„{pr)pdp 



(3) 



where F{p) is the Hankel transform oij{r) and we write 
the generalized-shifted function as /(?"l?"o) =R''°f{r), as a 
reminder that the shifted function is now a function of 
r and Yq. The operator 7?*^° acting on the function J{r) 
indicates a shift mfir) by Tq. The intuitive definition of 
the Hankel shift (generalized shift), as defined in (3), is 
that of the inverse Hankel transform of F(p)J„{pro), 
whereas the unshifted function would be the inverse 
Hankel transform of F(p) alone, without the multipli- 
cation by Jnipro). In essence, Equation (3) says that multi- 
plication by Jnipro) in the Hankel domain is equivalent to 
a generalized shift in the spatial domain. 

This definition follows the same expression used for 
the Fourier transform, that is: 



2n 



(4) 



In (4), the star denotes the complex conjugate and fol- 
lows the definition of generalized shift as given by Levi- 
tan. As previously pointed out, the simple shift, /(t- to), 
follows from the definition because e'"'e"'"'° = e"^''"'"' . 
For the Hankel transform with Bessel functions, no sim- 
ple equivalent expression exists, but the general struc- 
ture of the shift operation for the Fourier transform 
(left-hand side of Equation (4)) is the same. 

It is noted that using the definition of Hankel trans- 
form given in (1), the shifted function 7?''"/(r) = /(r[ro) 
can also be written as 



R'-'fir) =/(r|ro) = J I Jf{x)J„{px)xdx \j n{pro)J n{pr)pdp 



The interpretation of the generalized shift operator 
from Equation (5) allows the shift to be seen directly 
as an operation on the original untransformed func- 
tion. The definition as given in (3) allows for better 
physical interpretation of the definition (in particu- 
lar in comparison with the familiar Fourier trans- 
forms) and also permits the simple proofs to follow 
for the shift and convolution rules. 



Generalized shift rule 

In keeping with the standard shift-modulation rule of 
the Fourier transform, the definition of the Hankel 
generalized shift allows for a similar rule to be derived 
for the Hankel transform. The Hankel transform of 
the generalized-shifted function is given by application 
of Equation (1) to the generalized-shifted function in 
(3) as 

H„(/(rjro)) = J { J F{x)J„{xro)Jn{xr)xdx )j„{pr)rdr 

0 I 0 

= j F{x)J„{xrQ) j Jn{xr)J„{pr)rdr xdx 

0 0 




(7) 



where the order of integration has been reversed. 
Using the orthogonality (closure) of the Bessel func- 
tions (Watson 1995; Abramowitz & Stegun 1964) 



Jn{xr)Jn{pr)rdr = -d{x-p), 

X 



(8) 



it then follows that: 



ln(f{r\ro)) 



Six — 0] 
F{x)/„{xro) xdx 



Pip)Jn{pra) 



(9) 



f(x) J J„(pra)Jn(pr)J„(px)pdpxdx 



(5) 



In other words, a (generalized) shift in the spatial do- 
main is equivalent to multiplication by Jnipro) in the 
Hankel domain or, if J{r) Fip) then y(r|ro) <» f (^)/„ 
ipro). This follows the same rule as for the Fourier 
transform where the Fourier transform of f{r - ro) is 
given byf (cj)e"'"''». 



Baddour SpringerPlus 2014, 3:246 
http://www.springerplus.eom/content/3/1/246 



Page 3 of 6 



Modulation rule 

Since the Hankel transform is self-reciprocal, a "modu- 
lation" rule similar to that of the Fourier transform can 
easily be derived. The Hankel transform of Jni^Po)^/) is 
given by 



(10) 



«(/«(%)/('")) = j f{r)Jn{rpo)Jn{pr)rdr 
= F{p\p,) = RPoF{p) 



where the definition of the shifted transform follows 
by replacing r, Tq, p,f with p, po, r, F in the definition 
(3) so that 



RP«F{p) = F{p\p,) = j f{r)J„{rp,)j„{pr)rdr (11) 



In other words, if j{r) <» F[p) then it follows that /„ 

{rp,)f{r)^F{p\Po) = R'^m- 



Convolution rule 

With the definition of the generalized shift, we define 
the Hankel convolution of two functions as: 



l*Hf){r) = J giro)fir\ro)rodro (12) 

0 



0 0 ^, 0 



t-Uig*HfKr)) = J Jg{ro)f{r\ro)rodroI„{pr)rdr 

0 0 

= J Jg{ro)i^ jF{x)Jn{xro)J„{xr)xdx^ r^drQ} „{pr)rdr 

(13) 



Interchanging the order of integration gives: it then 
follows that 



^^n[[g*Hf)(r)) = j j j g{ro)J„(xrQ)rodro J„{xr)J„(pr)rdrF{x)xdx 



= j G{x)F{x) j Jn{xr)Jn{pr)rdr xdx = G{p)F{p) 
0 0 



--S[x-p) 



(14) 



In the preceding Equation, (14), the definition of the 
Hankel transform of g{r) has been used, in addition to 
the orthogonality of the Bessel functions. Equation (14) 
clearly states that the Hankel transform of the Hankel 
convolution is the product of the Hankel transforms, 
again in parallel with the standard result of Fourier 
transforms. Furthermore, interchanging g and / in the 
proof would give the same result, therefore it also fol- 
lows that (g* H f){r) = {f* H g){r) and that the Hankel 
convolution commutes. Therefore, we have that: 



{g*Hf){r) = {f*Hg){r)<>F{p)G{p) 



(15) 



Equation (12) is in keeping with the typical defin- 
ition of a convolution in the radial domain, where 
the simple shift f{r - vq) has been replaced with the 
generalized shift f{r\ro). Furthermore, we use the 
notation h to denote that this is a Hankel convo- 
lution, meaning that the generalized Hankel shift 
operator is used instead of the simple shift oper- 
ator. It is noted that other authors define a Hankel 
convolution without reference to the generalized 
shift operator. In all those cases, the integral of a 
triple product of Bessel functions is used to define 
the Hankel convolution, for example in (Tuan & 
Saigo 1995; Malgonde & Gaikawad 2001; de Sousa 
Pinto 1985; Belhadj & Betancor 2002). The math- 
ematical properties of Hankel convolutions are ana- 
lyzed in (Tuan & Saigo 1995; Malgonde & Gaikawad 
2001; de Sousa Pinto 1985; Belhadj & Betancor 2002; 
Betancor & Marrero 1993; Betancor & Marrero 1995; 
Cholewinski & Haimo 1966). 

We now proceed to find the Hankel transform of the 
Hankel convolution as defined in Equation (12): 



Multiplication rule 

Since the Hankel transform is self-reciprocating, then 
interchanging the roles of j{r) and g{r) with their trans- 
forms F{p) and G{p) in the previous derivation gives the 
result that the Hankel transform of the product /(r) g (r) is 
the Hankel convolution of their transforms {G'" h F)(p) = 
{F*hG){p), so that: 

f{r)g{r)^{G*„ F){p) = {F*„ G){p) (16) 
Example 

In this section, we apply the preceding definition to a 
commonly used function. The Boxcar (or gate) function 
is defined in Hankel frequency space as: 

The zeroth order inverse Hankel transform of the Boxcar 
function is given by: 




na{p)Jo{pr)pdp = I jQ{pr)pdp = -Ji{ar) (18) 



The function 2Ji{x)/x is often termed the jinc or som- 
brero function and is the polar coordinate analog of the 
sine function. Thus, we have that "Ji{ar)<^na{p) are a 
Hankel transform pair. Plots of the jinc and its Hankel 
transform box car are shown in Figure 1 and Figure 2. 

The generalized-shifted jinc function is thus given by 
Equation (3) and can be found in closed form via inte- 
grals given in (Watson 1995) as: 



-Ji{ar) 



I 



= / F{p)Jo{pb)Jo{pr)pdp = 



a 



(pb)Ja(pr)pdp 



a(rJa(ab)Ji(ar)-bJi(ab)Ja{ar)) 
b^-r'^ 



(19) 



A comparison of the original jinc function [a/r/i(ar)], 
its generalized shift i?''[fl/r/i(«r)] and its simple shift 
\^} \{a{r-b))\ for = 2 and ^ = 2 is given in Figure 3. 



Clearly from Figure 3 (and it should be obvious from 
Equation (18)), the generalized and simple shift are quite 
different and it would be a mistake to think that one can 
be used in place of another for modeling and simulation 
purposes. A comparison of the original jinc function, its 
generalized shift, and its simple shift for a = 1 and b=l/2 
is given in Figure 4. For smaller values of the shift, the 
generalized and simple shifts are closer to each other. 

Discussion 

The primary utility of the generalized shift function would 
appear to be that it is the function that permits the stand- 
ard shift, modulation, multiplication and convolution rules 
to apply when using the Hankel transform. 

We demonstrated in prior work (Baddour 2009; Baddour 
2011) that the Hankel transform is part of a two- 
dimensional Fourier transform in polar coordinates and 
that taking a convolution over the radial coordinate only 
(a convolution using the simple shift and only over r) does 
not lead to any simplification because the process of shift- 
ing over the radial coordinate destroys radial symmetry. In 
fact, it was shown in (Baddour 2009; Baddour 2011) that a 
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■ ■ Original fimaion gmgalized shift simple shift | 

Figure 3 Comparison of original function, generalized shift and simple shift, a = 2, b = 2. 



radially symmetric function shifted by 7 = (ro, 0o) in polar 
coordinates is given by: 

^{^-t) = Yl e*"-''»*y /(«) j Joipu) Jkipr^Vkipr) pdp udu 

(20) 

Thus, for a function shifted along the radial axis {do = 0), 
and supposing we are only interested in its values on the 
radial axis {6 = 0), then the simple shift in terms of the 
unshifted function is given exactly as: 



0 0 



(21) 



Generally speaking, the Hankel transform alone 
(without an accompanying angular coordinate Fourier 



transform to turn it into a two- dimensional Fourier 
transform) is most used in physical systems that have 
radial symmetry. Once the system is shifted in the radial 
direction - as is necessary to take a convolution - radial 
symmetry is lost and thus the proper, physically- 
meaningful transform would be a full 2D Fourier trans- 
form in polar coordinates. We showed in (Baddour 
2009; Baddour 2011) that if a full (radial and angular) 
shift is taken in defining the convolution, then the 2D 
Fourier transform in polar coordinates does possess the 
standard shift, modulation, multiplication and convolu- 
tion rules. 

What we have demonstrated in this paper is that if it 
is desirable to use only the Hankel transform and work 
with only the radial coordinate, then the Hankel trans- 
form does possess the standard shift, modulation, multi- 
plication and convolution rules - but only when the 
generalized definition of the shift is employed, not with 
the simple shift. 




10 15 

p^' • ■ Original function generalized shift simplt shift | 

Figure 4 Comparison of original function, generalized shift and simple shift, a = 1, b = 1/2. 
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Conclusions 

We have shown that the Hankel transform does possess the 
standard shift, modulation, multiplication and convolution 
rules - but only when the generalized version of the shift is 
employed, not with the simple shift. We demonstrated by 
way of a simple example that the generalized shift and sim- 
ple shift are not the same and thus not interchangeable for 
simulations of physical systems. The value of the general- 
ized shift is that it permits the standard Fourier rules to 
apply to the Hankel transform. For the purposes of calculat- 
ing physically meaningful convolutions, the simple shift in 
the radial coordinate should also be accompanied with 
an angular shift. A physically meaningful convolution im- 
plies integration with a shift over all physical coordinates 
allowed by the geometry of the problem. Fortunately, the 
full 2D Fourier transform in polar coordinates possesses 
the desired shift and convolution rules. 
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